On bulk viscosity and moduli decay 
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This pedagogically intended lecture, one of four under the header "Basics of thermal 

QCD", reviews an interesting relationship, originally pointed out by Bodeker, that exists 

between the bulk viscosity of Yang-Mills theory (of possible relevance to the hydrodynamics 

of heavy ion collision experiments) and the decay rate of scalar fields coupled very weakly 

to a heat bath (appearing in some particle physics inspired cosmological scenarios). This 

topic serves, furthermore, as a platform on which a number of generic thermal field theory 

("^ ' concepts are illustrated. The other three lectures (on the QCD equation of state and the 

^N) , rates of elastic as well as inelastic processes experienced by heavy quarks) are recapitulated 

►^ . in brief encyclopedic form. 

o 
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§1. Introduction and outline 

It has been a longstanding dream that experimental tests of thermal QCD 

through heavy ion collision experiments could yield theoretical insights that would be 

useful for some cosmological problems as well. These lectures covered selected topics 

Q^' within thermal QCD with this perspective in mind. The observables touched upon 



Oh: 



were the equation of state, shear and bulk viscosities, as well as the rates of elastic 
5^1 and inelastic reactions experienced by heavy quarks. Depending on the observable 

^^' the focus was either on elaborating on the basic concepts, on outlining the link 

between heavy ion collisions and cosmology, or on reviewing modern developments. 
^N I Because of severe page limitations, it is not possible to cover all of the lectures in 

(<— ^ ■ a detailed from in the proceedings. Rather, the choice has been made to concentrate 

Q> , on the second lecture, which is presented in some detail. The second lecture has been 

'^ I chosen because, first of all, the topic is rather elegant and concise, yet apparently 

. ■ little-known. Furthermore, it offers an opportunity to present a result which to my 

L^ \ knowledge is original, generalizing the analysis of ref.l^ from a scalar field theory to 

f~^ I a Yang-Mills theory. Finally, this topic can serve as a simple example permitting to 

introduce several tools (Kubo relations, Euclidean and Minkowskian correlators, the 
concept of transport coefficients) that are applicable much more widely in thermal 
field theory than to the very problem at hand. 

The presentation is organized as follows. In section [2] we review the current 
5^ \ theoretical understanding concerning the equation of state of QCD and the role that 

it may play in heavy ion collisions and in cosmology. Section [3l the main part of 
this presentation, concentrates on the topic announced in the title. The subject 
of section U] is the rate of kinetic thermalization of heavy quarks, caused by elastic 
scatterings with the light degrees of freedom in the hot plasma. Finally, in section [Sj 
recent developments in thermal heavy quarkonium physics, which to some extent 
can also be associated with the rate of inelastic processes felt by heavy quarks, are 
reviewed. 



typeset using PTPTeX.cIs (Ver.0.9) 



X 



2 M. Laine 

§2. Lecture 1: QCD equation of state 

In cosmology, it is in general an excellent approximation to set the baryon chem- 
ical potential to zero, so that thermodynamic potentials are functions of the temper- 
ature, r, only. Then the expansion rate of the system and the production rates of 
weakly interacting particles from the thermal medium (such as Dark Matter) depend 
on p(T),p'(T) and p"{T), where p{T) denotes the pressure^ The same functions 
(save without the contribution of photons and leptons, which have no time to ther- 
malize) dictate the expansion of the thermal fireball created in an energetic heavy 
ion collision. Hence the longstanding efforts, both analytic (via chiral perturbation 
theory at low temperatures, T <C 200 MeV, and the weak-coupling expansion at high 
temperatures, T ^ 200 MeV) and numerical (via lattice Monte Carlo simulations), 
to reliably determine these functions. At temperatures below one GeV, relevant for 
heavy ion collisions, the major qualitative finding, all but dashing early hopes of a 
spectacular scenario,'^'' has been that there is probably no phase transition at any T, 
only a smooth crossover.™ Nevertheless large-scale numerical efforts to determine 
p(T) and its derivatives for physical QCD in the infinite-volume and continuum limit 
go on. At temperatures above a few GeV, relevant for cosmology (only exceptional 
scenarios operate below this temperature ranged), the system should be address- 
able with weak-coupling techniques. Perturbation theory does suffer from serious 
infrared problems, however: besides proceeding in powers of al rather than asr*'' it 
also comes with non-perturbative coefficients starting at the order a^y^ That some 
coefficients are non-perturbative does not mean that they are beyond reach: indeed 
the non-perturbative inpulP can be determined via the use of effective field theory 
techniques.'^ As of today the non-perturbative part of the a^-term is known,'^ and 
it is possible to compile (systematically improvable) phenomenological values for 
p(T),p'{T) and p"{T) of the Standard Model, applicable in the whole temperature 
range from ~ GeV up to a few hundred GeV,*^ or even higher if necessary.!^ 

§3. Lecture 2: bulk and shear viscosities 

As a second example on the interplay between heavy ion collision experiments, 
cosmology, weak-coupling thermal field theory, and non-perturbative lattice simula- 
tions, we discuss observables known as the bulk and shear viscosities. 

3.1. Phenomenological background 

The bulk and shear viscosities are, in general, functions of the temperature and 
chemical potentials and parameterize gradient corrections to the energy-momentum 
tensor of a multiparticle system close to thermal equilibrium. At zeroth order in 
the gradient expansion one considers a system in which the temperature, T, and the 
four-velocity of the ("one-component") medium, u^, are constant in temporal and 
spatial coordinates; then the energy-momentum tensor, T'^^, is constructed out of 
symmetric and covariant structures proportional to the metric tensor, g^'^, and to 
u'^u'^ . At first order in gradients, the structure d'^u'^ + d'^u^^ can appear as well. 
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Thus, T'^'^ has the form (see, e.g., ref.'^) 

T"" = {p + ey^u"" -pg^'"' + Oii-n, Qid^u" + d^u")) . (3-1) 

Here p denotes the pressure, e the energy density, ij the "shear" viscosity, and C the 

"bulk" viscosity, and the metric convention (H ) was assumed. The equation of 

motion reads 9^T^^ = 0. 

More precisely, the shear viscosity coefficient rj is defined to be a function that 
multiplies the traceless part of the leading gradient correction; the bulk viscosity 
coefficient C multiplies the trace part. The explicit forms of the corresponding struc- 
tures are most simply displayed in a non-relativistic frame, \u^\ <C 1; then 

T^j ^{p-CV- v)6i.j - r]{div^ + djv' - ^6^jV • v) , (3-2) 

where V • v = diV^ and 5ij is the usual Kronecker symbol. 

Now, since the system generated in a heavy ion collision is relatively small, 
gradients can be large, with a scale given by the system size: 

( — ) ~ few fm ~ , (3-3) 

V " / ^QCD 

where Tqcd = 200 MeV. In the relativistic limit, rj ^ T^ andp ~ T^; this then implies 
that rjdu ^ pu^. In other words, gradient corrections could be as large as the leading 
terms. Setting aside the question of whether a gradient expansion can be justified 
at all in this situation, it seems conceivable that the viscosities could indeed affect 
the expansion of the thermal fireball created in a heavy ion collision. An example 
on the effects of the shear viscosity on the azimuthal anisotropics of particle spectra 
produced in a heavy ion collision can be found in refi--^ while possible implications 
of a (large) bulk viscosity have been explored in, e.g., ref.I^ 

We now turn to cosmology. An immediate qualitative difference between the 
matter created in a heavy ion collision and that filling the Early Universe is that the 
latter system is rather homogeneous. For instance, if we are interested in the overall 
expansion, then the system size can be identified with the horizon radius, which in 
the QCD epoch (T ~ Tqcd) is 

, "iPlanck 10^° 
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(3-4) 



This can be contrasted with eq. ()3-3p : thus gradient corrections are bound to be in- 
significant as far as the overall expansion is concerned. (Nevertheless, there may be 
other physics questions for which they do play a role; indeed some inhomogeneities 
do exist within the horizon, as is famously indicated by anisotropics in the cosmic 
microwave background and as is required for successful large-scale structure forma- 
tion. In spite of the fact that at early times their relative magnitude is of order 
~ 10~^ rather than of order unity, viscosities might still affect the evolution of den- 
sity perturbations; see, e.g., refs.'^''^^!) 

It turns out, however, that the bulk viscosity C makes a formal appearance in a 
completely different context.'^' In order to discuss this, we need to next recall some 
basic facts concerning the important role that scalar fields play in cosmology. 
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3.2. Scalar fields and the moduli problem in cosmology 

Although no fundamental scalar fields have been experimentally discovered up 
to date, they play a very important role in inflationary cosmology, and are also 
considered to be a generic feature of many models inspired by string theory or su- 
persymmetry. At an early stage, it is assumed that a significant amount of energy 
may be stored in the scalar fields. Later on, as the Universe enters the radiation 
dominated epoch, the scalar fields should "decay" , i.e. transform their energy to the 
known particles (electrons, photons, etc); this is called a reheating period. It is pos- 
sible to imagine scalar fields, however, which are coupled so weakly that they decay 
later than the inflaton field, or not at all. In the case of a delayed decay, one may 
use the scalar field to generate a possibly observable "second" component of density 
perturbations; such a scalar field is often called a "curvaton". Yet if the decay is too 
slow and the scalar field does not get rid of its energy density at all, it eventually 
becomes a dominant component, in conflict with observation. Such a situation is 
met particularly in connection with scalar fields called moduli, and the problem of 
their slow decay is the "moduli problem" . 

We note, in passing, that there are many other similar "dangerous relics" in 
cosmology: in particular theories predicting the generation of topological defects in 
the form of domain walls or monopoles are practically excluded, if the associated 
energy scale corresponds to physics beyond the Standard Model. 

To be concrete, let us consider a toy model in which a scalar field, (/?, couples to 
Yang-Mills fields: 

^ = l^i-a - ra')^ - -^F'^^'^F;, - ^ X F^^^F^, , (3-5) 

where m is some "small" mass scale, perhaps m ~ rn^vsY, while M is some large 
scale, perhaps M ~ ^n-pianck- The non-renormalizable coupling between ip and the 
Yang-Mills fields is suppressed by the heavy mass scale; hence the coupling between 
the two sets of fields is very weak. In ref.^^ a similar system was considered, however 
the scalar field if was coupled to another scalar field, denoted by x- 

If we now assume that the initial state of the system is such that ((/'(O)) 2> m, 
then the initial energy density in the scalar field is e^ ^ m^. According to eq. (j3-5|) 
the scalar field can decay to gauge bosons; however the coupling is suppressed by 
M, so the (vacuum) decay rate is i"' ~ m^/M^. The corresponding time scale, 
ty, ~ M'^/m^, can be contrasted with that in eq. (|3-4p : it is known from observation 
that a conventional radiation-dominated expansion is needed at least starting from 
the nucleosynthesis epoch (T ~ Tn = 0.1 MeV). (For completeness, a cartoon of a 
standard cosmological scenario is attached as fig. [TJ) If Af ~ w-pianck) the corre- 
sponding constraint t^ ^ t^ translates to m » ("T-pianckT^)^ ~ 10 TeV. While not 
impossible, this looks unattractive if m is to be related, for instance, to some lifted 
SUSY flat direction. 

3.3. Decay rate in thermal field theory 

The question now is, could the vacuum decay rate {F ~ tji'^/M'^) be modifled 
by thermal corrections associated with the "normal" degrees of freedom, represented 
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Fig. 1. A cartoon of a possible cosmological scenario. The age of the Universe, t, and the temper- 
ature of the plasma filling it, T, are related through t[s] ^ (T[MeV])~^ (cf. eq. 



by -F^jy? In particular, could the thermal plasma be "viscous" enough to dissipate 
the energy related to the oscillations of the "dangerous relic" cp well before the 
nucleosynthesis epoch, rendering ip harmless even if the mass m were only in the 
TeV range or below it? 

In order to answer this question, let us consider the equation of motion satisfied 
by ip. In the temperature range m <^T <^ M, the scalar field is "light" and oscillates 
slowly compared with the typical frequencies of the thermal degrees of freedom. In 
this situation, we can assume that the "fast" thermal fields see (p essentially as a 
constant, and have time to adjust to any given value. From their perspective, then, 
the system remains in thermal equilibrium, only in the presence of a background 
field ip. From the perspective of ip, on the other hand, the net effect of the thermal 
medium is to exert some friction which slows down its oscillations. Focussing on the 
infrared modes of <p (the modes with the slowest temporal and spatial variations), 
we can expect an equation of motion of the form 



Dv? + V^si^) = -Fip + OCif, VV, ^^ {y^f 



(3-6) 



The main task then is to determine the coefficient F (although we will also discuss 
thermal corrections to V^fj, particularly through a mass parameter w-gg). 

Now, a standard tool in thermal field theory is a "Kubo formula" , which indeed 
allows to determine "response" or "transport" coefficients of the same type as F. 
Note that F as defined by eq. ()3-6p is, by construction, a constant. It is given by 



(3-7) 



r=iim^(^ 



6 M. Laine 

where, for Ci^t = —^H\a\.-, /? is a spectral function related to the operator ^int: 

p{^,Q) = I e*"* (\ [^int(t,x),iJint(0,0)]\ . (3-8) 

The expectation value is taken with respect to the density matrix of the "normal" 
degrees of freedom which, as alluded to above, can be assumed to be thermalized. 



Sketch of a proof of eq. {3-7) 



Let us consider the form of eq. ()3-6p in Fourier space, with f oc e *'^*+*p'' c^; 

[-u^ + p2 + mis - i^r + 0(w^ wp2)]c^ = 0{ip^) . (3-9) 

We compare the dispersion relation from here with the "pole" appearing in the 
Euclidean propagator of (p after the analytic continuation a;„ — )• —i[uj + iO^): 

(3-10) 



-^n + P^ + ^E -w^ - «wO+ + p2 + Re 7l£; + i Im 7l£; ■ 
Taking p — )■ and denoting m^g = Re7l£;(— i[w + iO+], 0), we can identify 

r.-l,,'-'^^(-'l" + '°"l-«). (3.11) 

It remains to note that, as can be verified with basic path integral techniques 
in Euclidean spacetime, 

77sK, p) = - / e^-"-+*P- {H,^,{t, x) /7i„t(0, 0)) , (3-12) 

where < r < l/T, and that its imaginary part yields 

Im77£;(-z[a; + i0+],0) = -p(w,0) , (3-13) 

with p defined through the commutator in eq. (j3-8p : eq. ()3-13p will be justified 
presently. This completes the proof of eq. (|3-7p . 



Before proceeding to prove eq. (|3-13p . let us pause to contemplate an important 
point. The oscillation frequency given by eq. ()3-9p does not vanish, but is rather 
given by (for p — t- 0) 

M2 

where we have also included a thermal correction from the interaction in eq. 
Why is it then that F of eq. ()3-7p is evaluated at vanishing frequency? The reason is 
that p{uj,0)/lo shows generically a transport peak at origin, resembling a smoothed 
Dirac ^-function: it has the height F and a width, r/, determined by the scales of 



u;2«m,\~m2 + ^, (3-14) 
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the thermal system. For Yang-Mihs theory, parametricahy, ry ~ g^T}^ In the 
range m <^ T <^ M that we are interested in, however, rngfj ^ T according to 
eq. (|3-14|) . Therefore, if g is not exceedingly small, m^g <C ?? ~ g^T, and we can well 
approximate p{mcs,0)/m(,s by its limit Hui^^q p{u),0)/u), simplifying the problem 
and obtaining in any case the largest (i.e. most optimistic) F. 



Sketch of a proof of eq. i l3- J3|) 

Apart from the usual Heisenberg operators A{t) = e^^^A{0)e~'^^^ it is conve- 
nient to define "imaginary-time" Heisenberg operators as A{t) = e A{0)e~ , 
where < r < /3, /3 = 1/T. We furthermore define the correlators 

OD 

it e*^* 

oo 

OD 



n>{uj)= dte"^HA{t)A{Qi)), (3-15) 

J —oo 

n<{uj)= j dte*^*(i(0)i(t)) , (3-16) 

J —oo 



oo 
oo 



p{uj)^ / dte-*(iri(t),i(0) ) , (3-17) 



2 

oo \ ^ 



GE{r) = {A{t)A{^)^ , (3-18) 

fP . 
GE{uJn)= dTe"^"^GE{T); uj^ = 2TmT , n£Z. (3-19) 

Jo 

Here the expectation value is (...) = ^Tv{e^^^ {...)}, and it is easy to see 
that Ge{I3 — t) = Ge{t), i.e. Ge{t) is periodic. Spatial directions have been 
suppressed for simplicity but are trivial to incorporate. The correlator Ge{t) 
is theoretically a nice object, because it can be given a Euclidean path integral 
representation and thus a direct non-perturbative meaning. 

Inserting twice 1 = Yin l^)(^l ™ ^^^- P'15P - (|3-16p . it is straightforward 
to show that 11 ^{lo) = e~^'^7T>((x'). It then follows that 

piu;) = \ [n^iv) - 77<(u;)] = ^ (l - e'^^) n^{uj) , (3-20) 

or, conversely, that n^{uj) = 2p(a;)/(l — e~^'^). Making use of the fact that, 
through their definitions as correlators of Heisenberg operators, Gy{t) and 
Ge{t) are related by the analytic continution r o it, we finally obtain 



rp roo J 

GE{cOn)= dre'-"- / — e-*'^*77>(a;) 

t/ U L" — oo 

rl3 roo J 

Jo J-oo 2vr 



it^tT 







27r 1 - e-f^"^ 
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du p{u}) 


^{iui„—u})T 


TT 1 - e-^'^ 


iu>n — OJ 


do; p{uj) e"^'^ - 1 


TT 1 - e-^*^ iujn - OJ 


duj p{uj) 





(3-21) 

(To be precise, the validity of eq. ()3-2ip requires the presence of an ultraviolet 
regulator for spatial momenta, guaranteeing that both p{u}) and GEi^n) are 
cut off at large values of the argument.) Equation ()3-13p now follows by noting 
that the relation 

p(-] =fiTT6{x) (3-22) 



i^int = TjP'^^'Ku ■ (3-24) 



rc±^0^ 
implies that 

Im Ge {-i[oJ + iO+]) = p{uj) , (3-23) 

and by then inserting —UEiuJmP) from eq. (|3-12p for GE{^n)- 



3.4. Relation of the decay rate F to the bulk viscosity (^ 
Let us now consider the specific case of eq. ()3-5p , viz. 

1 

M' 

At this point it is good to realize that the structure appearing can be recognized as 
the "trace anomaly" of pure Yang-Mills theory, 

e = T^^^-^F-^^^F;,, (3-25) 

where 6o defines the /3-function related to the running coupling, 6o = llA'^c/3(47r)^. 

Sketch of a proof of eq. /i3-25\) 
With the convention 

S = Jdt |d3-2^x |-^F'^/^-f;,| (3-26) 

the energy-momentum tensor is 

^"^ = ^ (^5M.^'^"^^a/3 - F^%Ku^ ■ (3-27) 

In (5^^ = 4 — 2e dimensions its trace is 

^ 9b 



2-.2+-J7^77^^2+0(-) • (3-29) 
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The inverse of the bare gauge couphng {g"^ = 9^ — ^\f^T§;^ + C>{g^)) is 

1 1 ll^c 1 .^ (g^ 
gl g^^ 36 (4vr)2+''Ve 

So, in the Hmit e — )• 0, 

r^^-i^^F-Z^F^,. (3-30) 



Returning to eqs. (l3^2iH . ([3^ . we see that Fint = F^^'^F/^^/Af « -20/boM. 
Therefore eqs. (|3-7|) . (|3-8|) advise us to determine the transport coefficient related to 
the trace anomaly. However, as mentioned after eq. ()3-ip . the trace of the energy- 
momentum tensor is intimately related to bulk viscosity: in fact the latter can be 
expressed through another Kubo relation, 

C = ^lim- / e-*(M0(t,x),0(O,O)j\ . (3-31) 

Furthermore, motivated by the hydrodynamics of heavy ion collision experiments, 
the weak-coupling expression for C, has been worked out,'^^' 



hh'T^ 
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As the appearance of ln(l/as) in the denominator suggests, the computation is very 
non-trivial and necessitates a systematic resummation of the perturbative series. 
(Taking the Euclidean formulation as a starting point, there are also attempts at 
a non-perturbative determination of C in the temperature regime where the weak- 
coupling expansion can no longer be justified. '23) 

In conclusion, combining eqs. (|3-7p . ()3-8p . ()3-3ip . ()3-32p and the relation Hi^t ss 
—20/boM, the vacuum decay rate of the moduli fields, F ~ m^ / M"^, is overtaken at 
m ^ T ^ M by a thermal correction: 



r=lim- [ e'^'l- 
'^^0 w A,x \ 2 



^int(t,x),iJint(0,O) 



\2 rji3 ^3 



4x 9C (127ra^)^ T^ m 

^ l^M^ ^ ln(l/a,) M2 ^ M^ ' ^^'^^' 

In other words, the fact that there is already a plasma present "facilitates" the 
dissipation of the energy in the scalar field ip into that in normal radiation. 

3.5. Implications for the moduli problem 

We finally contemplate whether the rate in eq. p-33p could be fast enough to 
solve the cosmological moduli problem. Since the rate increases with T, we first 
need to figure out the initial temperature at which the oscillations start. To do this 
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properly we would need to rewrite eq. ()3-6p in an expanding background, but the 
upshot is that oscihations start once their frequency exceeds the Hubble rate (the 
inverse of eq. (|3-4|) ): 



mes ^ ^ T< Jrucs "ipianck ~ VmM . (3-34) 

"ipianck 

Then we should integrate the energy loss equation all the way down to the temper- 
ature T ~ m, at which point the vacuum decay takes over. We note, however, that 
within all of this range the decay rate falls far below the Hubble rate: 

(127ra,)2 r3 r2 

ln(l/a,) W ^ M ^^'^^> 



because, according to eq. (|3-34p . T/M < ^JrajM <C 1. Therefore, the decay is so 
slow that it essentially does not have time to take place within the lifetime of the 
Universe; thermal corrections quite probably cannot solve the moduli problem.^'' 

3.6. Summary 

The purpose of this lecture has been to illustrate various generic tools of thermal 
field theory, as well as the intriguing fact that systematic heavy ion collision inspired 
computations may find "exciting" applications in totally unexpected cosmological 
contexts. In the present example we could learn this way that the moduli problem 
remains a severe constraint even in the presence of thermal corrections, a fact to be 
taken into account in cosmological model building. 

§4. Lecture 3: heavy quark kinetic thermalization 

The production of a heavy quark and antiquark through gluon fusion, which 
then fly apart in opposite directions, is one of the most basic processes in a hadronic 
collision. In the heavy ion case, the "flying apart" part is non-trivial, however, due 
to the presence of a thermal medium which the heavy quarks have to surpass. The 
process is akin to Brownian motion, whereby the heavy quarks gradually lose their 
kinetic energy through collisions with the medium particles.^^''^^' This physics can 
be referred to as heavy quark kinetic thermalization, diffusion, momentum diffusion, 
drag, jet quenching, stopping, or energy loss. The situation becomes particularly 
tractable if we focus on a "late stage" of the process, in which the heavy quarks are 
non-relativistic with respect to the heat bath; then it can be described by Langevin 
dynamics, with the role of the stochastic noise being played by the colour-electric 
Lorentz force. Through a QCD generalization of the fluctuation-dissipation the- 
orem the problem thus boils down to the consideration of the 2-point temporal 
correlation function of the colour-electric field strength (rendered gauge-invariant by 
time-like Wilson lines) .l^'l^ The "transport coefficient" extracted from this corre- 
lator, conventionally referred to as k, has been the subject of some recent interest. 
A leading-order weak-coupling resuH!^ has been supplemented by a next-to-leading 
order correction,'2SJ which has however been shown to be so large as to question the 
validity of the weak-coupling expansion. Numerical simulations have been carried 
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out within so-called classical lattice gauge theory, confirming the existence of large 
infrared effects.'^ Computations through AdS/CFT techniques for similar processes 
in strongly coupled AA = 4 Super- Yang- Mills theory also suggest a much larger k than 
expected in leading-order QCD.^^^ '^^' All of this makes a strong case for attacking 
the problem with lattice simulations, a challenge that appears technically simpler 
than in the case of many other transport coefficients such as viscosities!^ If an 
answer can be obtained, it can be embedded in hydrodynamic simulations of heavy 
ion collisions (e.g. ref.'^SI') and eventually compared with experimental observations 
concerning the spectrum and azimuthal anisotropy of the decay products from heavy 
quarks.'Sy In cosmology, an analogous kinetic thermalization also plays a role, given 
that many Dark Matter candidates kinetically decouple in a non-relativistic regime, 
and their momentum distribution dictates the kind of structures that can form.'^ 

§5. Lecture 4: quarkonium in hot QCD 

Quarkonium physics in heavy ion collisions is somewhat similar to heavy quark 
physics discussed in the previous section, however in many ways more complicated: 
initial production is less likely because two heavy quarks and two heavy antiquarks 
need to be generated, and furthermore a quark and antiquark have to be in a suit- 
able kinematic range in order to bind together; in addition propagation through 
the medium is affected by several processes, such as scatterings experienced by the 
quark-antiquark colour dipole, decoherence of the quantum-mechanical bound state 
caused by interactions with the medium, as well as Debye screening of the potential 
that binds the system together. Recently, some progress has been made in sys- 
tematizing the study of these phenomena, by formulating the problem (implicitly 
or explicitly) in an effective field theory language.l^^ The setup makes use of scale 
separations, such as a^M ^ UgM ^ M at zero temperature, UgT <C al^^T <C T at 
finite temperature, and an additional UgM <C T <C al M to relate the two.l^ (It 
remains a challenge to promote the setup to the non-perturbative level.) An example 
of a generic feature emerging from the effective theory approach is that a concept of 
a static potential can be defined (as a "matching coefficient"), however it is static in 
Minkowskian rather than in Euclidean time. It is also complex unlike at zero tem- 
perature, with an imaginary part responsible for decoherence or, in frequency space, 
for the width of the quarkonium peak in the corresponding spectral function. Thus 
the static potential is different from that traditionally extracted from finite-T lattice 
QCD, which is purely real and relies on gauge fixing to the Coloumb gauge. Recently, 
an interesting attempt has been launched to extract the proper real-time potential 
from the lattice,'^ however further work is needed before conclusions can be drawn. 
Another insight is that purely perturbative studies may converge much better than in 
general (cf. the previous section),!^ because only short-distance physics plays a role 
for heavy quarkonium. Ultimately, one of the goals of these efforts is to determine the 
spectrum of the dilepton pairs produced from thermal quarkonium decays^--' (given 
that quarks disappear we may refer to this as an inelastic process). A cosmological 
analogue may exist in near-threshold two-particle Dark Matter production or decay, 
relevant for some of the most popular scenarios .'^''221 
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